Rotation of atoms in a two dimensional lattice with a harmonic trap 
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Rotation of atoms in a lattice is studied using a Hubbard model. It is found that the atoms are still 
contained in the trap even when the rotation frequency is larger than the trapping frequency. Fhis 
is very different from the behavior in continuum. Bragg scattering and coupling between angular 
and radial motion are believed to make this stability possible. In this regime, density depletion at 
the center of the trap can be developed for spin polarized fermions. 
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Recent spectacular progress in manipulating neutral 
atoms has opened the way to the simulations of com- 
plex quantum systems of condensed mater physics, such 
as high-Tc superconductors, by means of atomic sys- 
tems with perfectly controllable physical parameters. 
The systems studied include Bosc-Einstein condensations 
(BEC) for both atoms and molecules, paired states for 
fermions (BCS), the crossover between BEC and BCS 

So far, most of the experimental studies involving cold 
atoms were conducted in continuum. One of the major 
goals of studying ultracold gases in optical lattices is to 
understand the physics in condensed matter systems. In 
addition to topics already investigated in continuum sys- 
tems there are many effects of interest lately: Supcrfluid 
to Mott insulator transition [l3j , Bloch oscillation of par- 
ticles in lattice due to Bragg scattering 1A\, parametric 
atomic down conversion in BEC [l5j . and so on. 

In this letter, we will address the effect of a lattice on 
the rotation of atoms in a harmonic trap. Many phe- 
nomena of rotation atoms in a continuous trap have al- 
ready been investigated, both experimentally and theo- 
retically: the appearance of vortices [r| [Tt} in BEC and 
BCS sam ples PJJ, quantum Hall states for fast rotating 
fermions [la ], and vortex lattices in the lowest Landau 
level for BEC [HI, EH • On the other hand, lattices lead 
to many new effects under rotation, such as structural 
phase transitions of vortex matter [2(|. Also, near the 
superfluid-Mott insulator transition, the vortex core has 
a tendency toward the Mott insulating phase [2l|, and 
second-order quantum phase transitions between states 
of different symmetries were observed at discrete rota- 
tion frequencies [22| . 

In particular, it is well known that in a continuum the 
centrifugal force prevents atoms from rotating beyond the 
harmonic trap frequency That means that if the rota- 
tion is too fast the centrifugal force lets the atoms escape 
the trap. Therefore, a quadratic-plus-quartic potential 
was assumed to prevent the atoms fro m fly ing away from 
the trap at fast rotation frequencies [23|. In a lattice, 
however, as will be shown in this paper, it is possible for 
atoms to stay in the trap even if the rotation frequency is 
larger than the harmonic trapping frequency and density 
depletion at the center of the trap can then be developed 
for such a regime. 



For completeness, we first review the rotation of a par- 
ticle in a two dimensional (2D) continuum. In the rotat- 
ing frame, the Hamiltonian for the particle is 
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with L z = —id/dcf,, lu is the trapping frequency and f2 
the rotation frequency. Throughout this paper, units 
h = m = 1 are used, where m is the mass of the par- 
ticle. Eq. {l} is formally identical to the Hamiltonian of 
a particle of charge one placed in a uniform magnetic field 
2Qz and confined in a potential with a spring constant 
uj 2 — O 2 . This equation can be solved by separation of 
variables (r and </>) in polar coordinates where the radial 
equation is given as 



2p dp^ 'dp 



or 
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This decoupling of radial and angular degrees of freedom 
is different from the motion in a lattice, which is to be 
discussed later. For Q < lu, Eq. fl} has eigenvalues [13] 



E jtk = lu+(lu - fi)j + (lu + fl)k 



(3) 



where j, k are non-negative integers. Note that the angu- 
lar momentum states are eigenstates [13] such that there 
are only level crossings when f2 changes. At f2 = ui, the 
ground states become infinitely degenerate lowest Lan- 
dau levels. 

When f2 is bigger than lu (uj' 2 = Q 2 — uj 2 > 0), at 
r — > oo, the radial equation becomes -0" + r 2 Lo' 2 ip r = 0, 
which gives the non-vanishing asymptotic solution tp r cx 
exp [ztito'r 2 /2}. However, because the harmonic trap po- 
tential at r — > oo approaches infinity, the wavefunction at 
r — > oo has to vanish. This contradiction indicates that 
at SI > uj, Eq. (Q} has no solution. This means that the 
atoms would leave the trap because the centrifugal force 
2 r exceeds the restoring force —ui 2 r in the xy plane. 

Numerically, a system is treated by necessity as having 
a finite size. Thus a procedure needs to be determined 
to distinguish a "non-existing" wave function (meaning 
that the particles have left the trap) from a wave function 
where the particles are held in the trap. To this end, we 
discretize the simple analytic problem of Eq. ([T]) using a 
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FIG. 1: Energy as function of rotation frequency f2 in both 
the lattice and the discretized continuum. The energy in the 
lattice is shifted by At, the energy at the bottom of the band. 
In a real continuum, the atoms fly away from the trap at 
Q, > to (see text). 



central difference scheme with mesh size h. The results 
can be visualized like in Fig. [I] For well-contained parti- 
cles the analytic (Eq. ([3])) and numeric (Fig.Q} solutions 
agree in that the ground state energy E g does not depend 
on the rotation frequency f2. As soon as the particles are 
not contained anymore (since the centrifugal force is now 
stronger than the trap) the ground state energy plunges 
immediately. As can be seen in Fig. [I] this plunge hap- 
pens slower for larger mesh sizes. What is important here 
is that Eg depends in this (unphysical) case strongly on 
the mesh size. A similar case can be made by looking at 
the probability density |Vv| 2 of the particles in the trap, 
as shown in Fig. [2] For the case of a rotation frequency 
too large to contain the particles in the trap, it can be 
seen that only the fine-gridded solution (here with the 
mesh size h= 0.25) approaches our physical understand- 
ing of the particles being driven outward by putting the 
probability density all to the boundary of the numerically 
available space. Therefore, the wavefunctions depend on 
the mesh size only if the atoms are not contained and 
thus is another good measure (along with a varying E g ) 
for this situation. 

We now study the case of zero-temperature atoms in 
a lattice placed inside a harmonic trap. We will see that 
in this case, the particles stay contained even when the 
rotation frequency fl is slightly larger than the trapping 
frequency uj. In the rotating frame, the single band Hub- 
bard Hamiltonian is 1221 
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(i,3) 



-E^)"* ( 4 ) 



where indicates a sum over nearest neighbors, 

V(fi) = £l 2 rf/2 is the harmonic trap potential, t is the 
hopping term that describes tunneling between neighbor- 
ing sites, rii = cfcj is the number operator with cf (c.i) 
the fermion creation (annihilation) operator at site i, and 
H.c. means Hcrmitian conjugate. In this paper, t is set 



FIG. 2: Diagonal density profile as function of the radius r 
from the center of the trap for the discretized continuum. 
When the rotation frequency Q is 0.11, which is just above 
the trapping frequency u = 0.1, the finer the mesh size h, the 
closer the wave function are pushed towards the boundary by 
the centrifugal force. The inset, accordingly, shows that the 
energy drops when the mesh gets finer. 



to one and used as a reference unit for uj and tt. flK rep- 
resents the centrifugal term with JQj = r^r^ sin a^j/d 2 , 
where denotes the distance from the axis of rotation to 
the ith site, onj is the angle subtended by the zth and jth 
sites with respect to the axis of rotation, d is the lattice 
constant, and (3 = 0.493 is the dimensionless constant 
characterizing the lattice geometry and depth [22|, H3|- 
We note that, different from that in continuum, the ef- 
fective mass in the lattice n(k) = (2td 2 cosfcd) -1 depends 
on momentum k and may become negative. 

To avoid confusion here it should be noted that with- 
out a trap, the particles will immediately fly away in- 
dependent of the magnitude of f2. This means that the 
lattice potential (in particular, the tunneling t term) has 
absolutely no confining effect on the atoms. As a result, 
the addition of the lattice potential to the trap potential 
cannot explain the above-threshold rotation which will 
be explained here. 

Figure Q] shows that in the presence of the lattice E g 
does not depend on the rotation frequency for Q < uj as 
it is the case for a continuum. For comparison, in Fig.[T] 
we choose the lattice constant such that n(k = 0) = 
1 and find that the O-dependence for Q, > uj of E g is 
the same as for the case of the "discretized continuum" 
(with the same effective mass [i = 1) by setting the mesh 
size h = 1/ y/2. The similarity in the energy dependence 
for the lattice and discretized continuum let us conclude 
that it is just a consequence of the discrete nature of 
the lattice. Besides, it is to be remembered that the 
infinite degeneracy of the ground state at il = uj in the 
continuum is lifted in the lattice and thus the physics 
associated with the lowest Landau levels will be different 
from that in the continuum. A detailed study of this will 
follow in the future. 

Figure [3] shows the density distribution at various 
51 around uj = 0.1. Two features deserve attention: 
First, the atoms still are contained at fi = 0.11 > 
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FIG. 3: Ground state density distributions at = 0.09 (a), 
0.105 (b) and 0.11 (c). As the rotation frequency SI becomes 
larger than the trapping frequency lu = 0.1, the ground state 
wave functions gradually change into the 4-lobe structure and 
the atoms are kept hold in the trap. For better visibility, the 
maximum density are normalized to one for each plot. 

uj because if we increase the lattice size, both the 
ground state wave function \ip^ and energies E g do not 
change. At larger rotation frequency SI, different lattice 
(SIZE/CONSTANTS) lead to different density distribu- 
tions and ground state energies, which we interpret as 
particles escaping the trap. In addition, the density at 
the center of the trap becomes depleted, which is im- 
possible for the ground state in the continuum as shown 
in Ref. [Uj]. When the lowest four states become de- 
generate, the ground state develops a 4-lobe structure 
(Fig. [3t). Moreover, when fl increases beyond ui, the 
ground states continue to change smoothly and move 
away from the center as shown in Fig. [3] the particles 
in the ground states only move away from the center to a 
distance determined by the rotation frequency, but don't 
leave the trap altogether. It is also found that with the 
same lu, the particles in the ground state of a smaller lat- 
tice leave the trap at smaller values of SI, as opposed to 
the universal threshold lu that is present in continuum. 
This difference between the lattice and the continuum is 
one of the major results of this paper. Note that the 
symmetric 4-lobe structure is determined by the primi- 
tive cell geometry's fourfold symmetry, not by the over- 
all lattice shape, because, for example, 150X100 lattice 
gives the same diagonal density distribution as 150X150 
at lu = 0.1 and SI = 0.11. In other words, the fourfold de- 
generacy reflects the discrete rotational symmetry of the 
underlying square plaquette in the lattice [22 |. Further- 
more, when the 4-lobe structure is finally formed, both 
at x = and at y = the densities approach zero. 

The atoms rotating in the lattice are always in mo- 
tion, so there are always currents between neighbor- 



ing sites, which are calculated using 
it{a,iOj — H.c.) + tt^id^ + H.c.) [22|]. One example of 
the current pattern is shown in Fig. 2] While the proper- 
ties of the 4-lobe states will be for future publication, the 
important thing here is that the persistent motion allows 
Bragg scattering to explain the stabilization of atoms in 
the lattice. 
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FIG. 4: The current distribution at one of the lobe in Fig. |3j;. 
A center for the current pattern is clearly seen at site (-10,-10) 
by nearly disappearing current at that site. However, Fig. [3fi 
shows that the density at this site is nonzero, different from 
a vortex core. 



To appreciate the role of Bragg scattering, we recall 
the Bloch oscillations in the lattice [3] : the driving volt- 
age across the lattice does not produce a net current for 
the electrons, instead it produces periodic current, as the 
Bragg scattering completely reflects the electrons. So it 
may not be so surprising to see that the Bragg scattering 
can hold the atoms in the trap even when S7 is slightly 
larger than lu. Furthermore, as the continuous symmetry 
is broken in the lattice, the correlated motion in radial 
and angular direction allows Bragg scattering between 
these directions. This is believed to also be the reason for 
the containment of the particles for rotation frequencies 
beyond SI = oj. Additionally, we note that Bragg scatter- 
ing is not limited to the Hubbard Hamiltonian Eq. ([4|, 
because the discretized continuum Hamiltonian Eq. (fT]) 
with large mesh sizes gives similar results (the dashed 
line in Fig. [1]) as that in lattice Eq. ((4]) . 

So far, we have shown that the quantum statistics of 
atoms does not play a role for the stabilization. Next, we 
point out that density depiction at the center of the trap 
for spin polarized fermions is possible. To this end, spin 
polarized fermions without interaction are considered in 
the lattice at zero temperature. The overall density of 
putting N rotating fermions into a lattice, according to 
Pauli's exclusion principle, is the summation of the den- 
sity of the lowest N states. Note that the very low ly- 
ing states all have zero population at the center of the 
trap, which causes the non-interacting fermions not only 
to have a density plateau at SI < lu (the solid curve in 
Fig. [5]) like in the continuum [l8| , but one also sees den- 
sity depletion at SI > lu as shown by dashed line and 
dot-dashed line in Fig. fSJ This is another major result of 
this paper. In comparison, the density depletion cannot 
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be realized in a continuum with a harmonic trap because 
the atoms fly away at Q > u. 



0.03 

C 0.025 
CD 

Q 0.02 









/ \ o-n nQ 




f \ £2=0.11 




/ \ £2=0.1 3 







•60 -40 -20 20 40 60 
X 

FIG. 5: Density distribution of 50 noninteracting spinless 
fermions at the cross section with y = 0. When £1 is close 
to but smaller than u = 0.1, the plateaus are developed re- 
flecting the underlying Landau-level wave functions [la ]. On 
the other hand, when f2 is larger than u, density at the center 
is depleted. 

In the above, coupling between the angular motion and 
radial motion as well as Bragg scattering were used to 
explain the stabilization. One may also think that the 
effective mass, which depends on k and could be nega- 
tive in the lattice, may also stabilize the atoms. However, 
the wave packets with fi(k) > will not be contained in 
the trap at fl > u> if only the effective mass contributes 



to the stabilization. Therefore, while it may play an im- 
portant role, the effective mass alone cannot explain the 
stabilization. 

In lattices, when ui is comparable to or greater than t, 
the coupling between the angular motion and radial mo- 
tion is enhanced. As this enhancement makes the Bragg 
scattering between radial direction to angular direction 
larger, the atoms will stay in the lattice with larger ro- 
tation frequency. This means that to experimentally 
demonstrate the above threshold rotation, larger trap- 
ping frequency is desirable. However, deep harmonic trap 
makes the single band Hubbard model not be reason- 
able for describing an optical lattice system [25[. There- 
fore, some intermediate trapping frequency is preferred 
for demonstrating the above-threshold rotation experi- 
mentally 

To conclude, the rotation of atoms in an optical lattice 
is studied using a Hubbard model. It is found that the 
atoms are contained in the trap even when the rotation 
frequency exceeds the trapping frequency, which is very 
different from the continuum case. Bragg scattering and 
the coupling of angular and radial motion make the above 
stability possible. In this regime, density depletion at the 
center of the trap can be developed for spin polarized 
fermions. 
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